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We have performed HNC calculations for dense beryllium plasma as studied experimentally using x-ray Thom- 
son scattering, recently. We treated non-equilibrium situations with different electron and ion temperatures 
which are relevant in pump-probe experiments on ultra-short time scales. To consider quantum effects ade- 
quately, we used effective pair potentials to describe the interactions. Results are compared with classical as 
well as quantum corrected Debye model calculations. 
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1 Introduction 

The dynamic structure factor S (k,oj) is the spectral function of the density fluctuations. It determines funda- 
mental properties of the plasma such as pair distribution functions, equation of state data, and transport coeffi- 
cients [ 1 1. Simultaneously, it gives direct access to the frequency spectrum of amplitude-modulated electromag- 
netic waves scattered off the plasma and, therefore, is a key quantity for plasma diagnostics. We are especially 
interested in strongly coupled plasmas with densities typical of solid state and temperatures of several eV which 
are also known as warm dense matter (WDM). Intense x-ray radiation sources are needed to perform scattering 
experiments in WDM targets. Free electron laser facilities like FLASH at DESY in Hamburg [2] will be available 
to deliver such radiation. Alternatively, pump-probe experiments are also planned at GSI Darmstadt using the 
future FAIR facility 0. 

Pioneering pump-probe experiments on warm dense beryllium have been performed at the Laboratory for 
Laser Energetics with the Omega laser in Rochester. They were successful in spectrally resolving the non- 
collective (particle) scattering characteristics of beryllium [4] and carbon [5]. Recently, it has been shown J6| 
that important plasma parameters such as electron density and temperature can be inferred from the collective 
scattering feature (plasmons) of the measured dynamic structure factor. 

Therefore, a detailed analysis of the dynamic structure factor S (k,u>) in WDM states is needed in order to 
account for strong correlations and scattering processes between the particles which determine the damping of 
excitations processes Hill- Time resolved measurements on ultra-short time scales would offer the possibility to 
study non-equilibrium states, e.g. different electron and ion temperature and the relaxation into thermodynamic 
equilibrium, see ll9l [lOl fTTI [ill . 

We present an approach describing a plasma consisting of electrons and ions having different temperatures. 
In order to calculate the static structure factor S (k) for dense beryllium plasma we solve an Ornstein-Zernike- 
like equation using a hypernetted chain (HNC) closure relation. Effective pair potentials are used to account for 
quantum effects. Results are compared with classical as well as quantum corrected Debye model calculations of 
Gregori et al. |9|. Especially, non-equilibrium situations with different electron and ion temperatures which are 
relevant on short time scales are studied motivated by the recent experiment |6). 
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2 Theory 

The rigorous way to describe a multi-temperature plasma is to start from the BBGKY hierarchy jT3l[T4l[T5l[T6l . 
We are interested in the two-particle distribution function which allows, in principle, to calculate most of the 
physical quantities, e.g. pressure or internal energy. Here, we are especially interested in the structure factor, in 
order to interprete the scattering results. 

We consider a homogenous and isotropic system where binary interactions V ab and the two-particle distribu- 
tion function F ab depend on the distance ri 2 = ri — r 2 of the two particles only. In absence of external forces 
the equation of motion for the two-particle distribution function F a b reads 

/ d dV ab (r 12 ) d dV ab (r 12 ) d pi 9 p 2 d \ 

Nt + a T^— + a * a -z—) Fab(ri2,Pi,P2,t) (1) 

\ot dri dpi ar 2 op 2 m a dri m b dr 2 J 

Efj3 j3 dV ac {r 13) d 
n c J d J r 3 d J p 3 — g^ F abc (ri,r 2 ,r3,pi,p 2 ,P3,t) 

Et ,3 j3 dV bc (nj3) d 
n c J d J r 3 cPp 3 — — F abc (r u r 2 , r 3 , pi, p 2 , p 3 , t) , 

where indices a, b, and c denote the plasma species. The quantities n c are the particle densities, and F a b c 
is the three-particle distribution function, which itself couples to higher-order functions. One has to find an 
approximation to truncate equation (Q]l. 

As mentioned before we are interested in a model of the plasma, where each species c is described within a lo- 
cal equilibrium with temperature T c . This is possible if the momentum relaxation for particles of the same species 
is faster than for particles of different species. This fact is, for example, fulfilled having mass ratios far away from 
one between different species. Then the momentum distribution separates into a product of independent Maxwell 
distributions with temperature T c . For the two-particle distribution function F a b one gets 

Fab (ri 2 ,pi,p 2 ) = g ab (1-12) F a (pi) F b (p 2 ) (2) 

with F a (pi) = exp (—0 a p1/ (2m a )) / (2nm a / f3 a ) 3 ^ 2 , where m c is the mass of the species c and [3 C = (ksT c ) 1 
the inverse temperature. Higher-order distribution functions are defined likewise. Using this ansatz in equation (Q]) 
yields 

a Pi dV ab (r i2 ) p 2 dV ab (n 2 ) pi d p 2 d \ 

Pa 5 h Pb 5 1 q h — ^— 5afc (ri 2 ) F a (pi) F fc (p 2 ) (3) 

m a ot\ mb oy 2 m a or\ nib or 2 I 



= -p a — >^«c / d r 3 g abc (ri, r 2 ,r 3 ) F a (pijFb (p 2 ) 

- Pb — ) ™c I d J r 3 ffahc(ri,r 2 ,r 3 )F a pi)F fc (p 2 . 

mb ^ J or 2 



Consequently, we concentrate on the two-particle spatial distribution function g a b only. We have first, in 
equation (O, to integrate over the momenta. Therefore equation (O is multiplied by pi/m a and by p 2 /mb, 
respectively, yielding two equations. Performing the integration over the momenta for each equation results in 

1 dV ab {r 12 ) , k B T a d 



dr 12 m a dr 



12 



9ab (r i2 ) (4a) 

> n c / d r 3 3af>c(ri,r 2 ,r 3 ) 

m a J dri 

1 dV ab (r 12 ) , k B T b 8 \ 

9ab{ri2) (4b) 

f ,3 dV bc {r 23 ) 

J r3 — dr~ 2 — 9abc ^ ri ' r2 ' r3 ^ ■ 



m b dr i2 m b dv X2 
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Summing up equations d4ab and ( Rbl leads to 

dhgb (rig) gKb (rig) , , 

^^ +/3afc ^^ 3afc(ri2) (5) 

fl "lb /"^3 ^Kc(ri3) , . 

= -Pab ; > ™c / a r 3 3 afcc (ri,r 2 ,r 3 ) 

m a + m b ^ J dvi 

+ Pab ; > « c / a r 3 g abc (n, r 2 , r 3 ) 

ma+mj^ 7 dr 2 

c 

with /3 a f, = (ksTah)^ 1 being the mass-weighted inverse temperature of two species a and 6 with T ab = (m a Tb + 
mbT a ) I (m„ + m&) and h a b = <? a b — 1 the total correlation function. 

Result © represents the equation of motion for the two-particle spatial distribution function g a b- As noted 
earlier it is coupled to higher-order contributions in the hierarchy via the three-particle spatial distribution function 
g a bc- In order to solve equation dD a truncation is needed. 

Assuming thermal equilibrium at the temperature T equation (O reads 

dKb (rig) dVgb (rig) 

— « hp — ^ gab(r 12 ) (6) 

OT12 OTi2 

TO6 



/"j3 fl ^Kc(ri3) , s 

2_^n c J d°r 3 (3 — g abc (n,r 2 ,r 3 ) 

m ° /"^3 «^6c(r 23 ) , s 

j > n c / d r 3 ^ g abc (ri,r 2 ,r 3 ). 

^a + m b ^ J dr 2 



It is known that equation © can be decoupled formally by introduction of the direct correlation function c ab , 
giving 

dh ab (r 12 ) dc ab {r 12 ) m b ^ f 3 dc ac (r 13 ) 

— 5 5 = 1 " c / d r 3 5 /i c & (r 23 (7) 

OT i2 dr i2 m a + m fc J Sri 

1 ) "c d J r 3 h ca (n 3 ) , 

m a + m b J dv 2 

which is equivalent to the Ornstein-Zernike equation IfTTl 

hab (ri 2 ) - c ab (r i2 ) = X! " c / rf3r3 Cac ( ri3 ) ^ ch ( r23 ) ■ 

To lowest order in the interaction, one has c ab ss h ab w —f3V a b- 

For the multi-temperature regime, it is reasonable that decoupling the hierarchy by a direct correlation function 
c a b fulfills, to lowest order, the relation c ab ps h ab ps —(3 a bV a b- Thus the corresponding multi-temperature 
Ornstein-Zernike equation reads 

dhab(ri 2 ) dc ab {r 12 ) m b ^ f 3 1 dc QC (r 13 ) 

= +/?ab ; > n c / d r 3 h cb (r 23 ) (9) 



dri2 dr 12 m a + m b ^ J f3 ac dri 

c 

a m a m-^ f 3 1 dc bc (r 23 ) 

- Pab ; > n c / tf 3 ^ ft ca (r 23 ) . 

m a + m b ^ J (3 bc dr 2 

In order to solve equation (O a closure relation is needed, connecting the direct correlation function c ab to the 
total one h ab . In equilibrium there exist various closure relations lfT8l[T9l [20 21 , 22 1 that can be easily transferred 
to a multi-temperature plasma by the same way described above, i.e. introducing the mass-weighted temperature 
Tab for the species a and b. We use a HNC-like closure relation written as 

Cab (ri 2 ) = ft-afc (ri 2 ) - lng a b (r i2 ) - pafcKb (ri 2 ) • (10) 
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The non-linear set of equations (O and ( fTob is solved in the usual way ll23l giving results for the direct and 
total correlation functions, c a b and hat,. The structure factor S a b is defined by the well-known relation 



Sab (k) = 8 ab + y/n a n b J d 3 ri 2 exp (-ikr 12 ) (g a b (V12) - 1) 
= Sab + y/n a n b Kb (k) . 



(11) 



3 Results and Discussion 

We have considered a two-component beryllium plasma consisting of electrons and ions with an effective charge 
state Z. For the description of the interactions V a b there exist numerous effective pair potentials, e.g. see 11241 |25l 
|26]|27]. For all these potentials, the main goal is to incorporate short-range interactions to avoid divergencies 
which result from the behaviour of the Coulomb potential V a b (712) = (qaqb) / i^Sor^) at zero distance. In 
classical schemes, this problem may be solved by hard-sphere-like potentials. An improvement is given by using 
pseudopotential theory as known from solid state physics. In a first simple step we apply quantum mechanics by 
evaluating the two-particle Slater sum S a b- One gets an effective potential V a b, see Morita ]28 1, via 



ib (ri2) = exp (-PabV a b (ruj))- 



The two-particle slater sum itself is defined as 

Kb (ri 2 ) = 8^\ 3 ab 



(-lf Sa S ab 



2s a + 1 



exp (-0abE a ) |* a (ru, 0, 0)| 2 . 



(12) 



(13) 



Here the thermal wavelengths X a b = fr\J Kb/ (2m a b) with m ab x = 



1 and the spin quantum number 



s a are introduced. The quantity a stands for the quantum numbers n,l, and m for bound states and k,l, and m for 
scattering states with fy a being the corresponding wave function, respectively. 

The first-order Born approximation at low degeneracy in equation [\3[ leads to the Kelbg potential 11241 



Vab (m) 



qaqb 



1 - exp [ --T2 - ) + V^- — erfc [ - — 
Kb 



ri2 r / ri 2 
7r- — eric I — 

Xab V Kb 



(14) 



with erfc (x) being the standard complementary error function. This type of potential has already been considered 
in El. 

For the electron-ion interaction one has to exclude the effect of bound states. We applied in the subsequent 
calculations a potential derived by Klimontovich and Kraeft [29 1 (a simplified version of their equation (42)) 



V ei (ria) 



16 



k B T et £,j 

i-rreo 



(15) 



where £ e j = (q e qifl e i) / {^£oKi)- In this potential, only the ground state of a bound electron is excluded. 
Both the Kelbg potential (TBI) and the Klimontovich-Kraeft potential ( TT3T > are finite at zero distance. The equa- 
tions (V2\ . ( Tol l, O, and $15[ , derived for thermal equilibrium, are generalised to a two-temperature situation by 
replacing the temperature T with the mass- weighted temperature T a b as explained in section [2] The interaction 
potentials V a h are displayed in figure Q] 

We calculated the ionic structure factor Su within the formalism derived in section[2] For the description of the 
electron-electron and ion-ion interaction we used the Kelbg potential (TT4b . Because we found the Kelbg potential 
not applicable for the electron-ion interaction in that low temperature region, we applied the Klimontovich-Kraeft 
potential (1151 , Results for the ionic structur factor Sa are shown in figure [2] versus the wave number k for a 
beryllium plasma with charge state Z — 2.5, ion density n,; = 1.21 x 10 23 cm~ 3 , and fixed electron temperature 
T e = 12 eV. We compare calculations for different ratios of electron T e and ion Tj temperature. Also shown are 
results by Gregori et al. [9| and for the Debye approximation for a pure ionic system, Sa (k) — k 2 / (fc 2 + n% ) 
with n 2 ^ = (Z 2 e 2 nA / (eoksTi), which can be considered as the limit for weak correlations. 
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Fig. 1 Effective pair potentials \ V a b \ used in the calculations for a beryllium plasma shown against the distance 
ri2 for charge state Z = 2.5 and fixed electron and ion temperature at 12 eV. Notations are as follows: solid 
line for the Coulomb potential, dashed line for the Kelbg potential dl4t . and dotted line for the Klimontovich- 
Kraeft potential d 1 5 b - The interactions are marked by a filled circle for the electron-electron, diamond for the 
electron-ion, and triangle for the ion-ion interaction. For the ion-ion interaction the Kelbg potential is close to 
the Coulomb potential, but has a finite value at zero distance. 




Fig. 2 Ion-ion structure factor Sa for a beryllium plasma with charge state Z = 2.5, ion density rn = 
1.21 x 10 23 cm -3 , and fixed electron temperature T e = 12 eV plotted versus the wave number k. The result of 
the present work is displayed as solid line. The work of Gregori et al. (9) and the Debye result are shown with 
dashed and dotted line, respectively. Three different temperature ratios T e /Ti between electron T e and ion T% 
temperature were examined marked as filled circle for T e /T,: = 1, diamond for T e /Ti = 2, and triangle for 
T e /Ti = 4. The Debye approximation is plotted for Ti = 12 eV. 

In contrast to the present paper, Gregori et al. [9] did not use the HNC equations, but generalised the screened 
pseudopotential approach of Arkhipov and Davletov [30] for a two-temperature system. In this approximation 
the ionic structure factor is given by Sa (k) = 1 — n^u (k) / (fcsT/). The effective ion temperature T[ = 
y/T? + 0.152 with Tjj being the Debye temperature accounts for ion degeneracy at low temperatures. The 
pseudopotential is a modified Kelbg-Deutsch potential Il25ll26l with Debye-like screening corrections. 

Consider in figure |2] the curves for the temperature ratio T e /Ti — 1: For small k, the present result is close 
to the Debye result. Our calculations show a small but nevertheless non zero value at k = 0. With increasing 
k, the ionic structure factor is increasing rapidly showing a peak at k w 1.8 Od . For values of k greater than 
k « 4.5 8j , the function is almost one, that means the system is uncorrelated at that scale. The result of 
Gregori et al. clearly shows a different behaviour. Their curve starts with a finite value around 0.3 at k = 
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indicating long range polarisation effects. With increasing k the function increases slower showing no peak and 
turning into a Debye-like behaviour. 

Concerning the temperature dependence there is a common trend in the results of the present work and those 
of Gregori et al. (5J. With decreasing ion temperature Tj, the value of the ionic structure factor Su for small k is 
lowered. This can be attributed to a weakening of the correlations between the electrons and ions. The system is 
approaching towards a one-component plasma (OCP) with static ions in a screening background. Within the HNC 
calculations this leads to a sharper and stronger peak. The occurrence of such structures in our calculations is a 
progress in the description of multi-temperature plasma resulting from the inclusion of higher-order correlations, 
see for comparison l23l . 

0.40 T 
0.30-\ 




0.00 J 1 1 1 1 1 1 1 1 1 1 

1 2 3 4 5 6 

T /T. 



Fig. 3 Ion-ion structure factor Su for a beryllium plasma as a function of the ion temperature Tj, The electron 
temperature T e is fixed at 12 eV, the ion density is rn — 1.21 x 10 23 cm -3 , charge state is Z — 2.5, and 
k — 0.54 . The result of the present work is displayed as solid line. The work of Gregori et al. (9| and the 
Debye result are shown with dashed and dotted lines, respectively. 

For the specific experiment [6 |, the ionic structure factor at k = 0.54 a^ 1 is of interest. For Thomson scattering 
at a wavelength of A = 0.42 nm, the k value corresponds to a scattering angle of 9 — 40°, i.e. forward scattering. 
For this k value the ionic structure factor Su is considered as function of the ion temperature Tj and of the charge 
state Z in figures [3] and |4] respectively. Figure [3] shows, for this small k value, a systematic decrease of the 
structure factor for decreasing ion temperature. The absolute values however are model dependent. 

The dependence on the charge state is plotted in figure |4] for T e = T, ■ = 12 eV. The curves look similar to 
figure |3] This can be understood easily keeping in mind that a lower ion temperature as well as an increased 
charge number Z are connected with a stronger coupling of the ions in the system. 

4 Conclusions 

We have calculated the static structure factor S (k) for solid-density beryllium consisting of electrons and ions. 
Effective pair potentials were used to include quantum mechanical corrections. We studied non-equilibrium 
situations for a two-temperature plasma. In this case a mass weighted temperature occurs. We find a strong 
influence of correlations on the static structure factor compared with simpler approaches such as the Debye 
model with quantum mechanical corrections of Gregori et al. J51 and the classical Debye model. Therefore, the 
correct treatment of all correlations in a two-component plasma is crucial for the description of non-equilibrium 
states and their relaxation into equilibrium. A consistent derivation especially of the effective electron-ion and 
ion-ion potentials, is, of course, needed and remains the subject of further work, e.g. accounting for finite size 
corrections. 
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Fig. 4 Ion-ion structure factor Su for a beryllium plasma as a function of the ion charge state Z. The electron 
temperature T e as well as the ion temperature Ti are 12 eV, the ion density is m = 1.21 x 10 23 cm -3 , and 
k = 0.54 ag 1 . The result of the present work is displayed as solid line. The work of Gregori et al. (5| and the 
Debye result are shown with dashed and dotted lines, respectively. 
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